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1 . INTRODUCTION 


The present practice of using high-bypass turbojet engines has 
resulted in a decrease in jet noise. However, these engines emit noise 
from their inlet nacelles above a desirable level. The present work is 
part of a considerable effort being made to reduce inlet noise. 

One promising approach to the reduction of inlet noise is the use 
of a high subsonic Mach number inlet, or partially choked inlet, in 
conjunction with an acoustic duct liner. The use of choked inlets has 
long been recognized as an effective means of reducing upstream propaga- 
tion although such inlets require careful design to prevent excessive 
losses in compressor performance. However, the physical mechanisms 
responsible for the noise reduction in high-subsonic Mach number inlets 
are not completely understood, and techniques for the theoretical 
analysis of sound propagation through regions of near-sonic mean flow 
are not available. Two major problems must be overcome in the develop- 
ment of such a model: (1) the mathematical techniques for the calcula- 
tion of sound propagation in ducts are wel 1 -developed for parallel ducts 
but are not fully developed for ducts of varying cross section that 
carry mean flows with strong axial and transverse gradients; (2) linear 

acoustic equations are inadequate to describe acoustic propagation in 

1 2 

regions of near-sonic mean flows ’ . In the investigation presented 
here, the first of these two problems was addressed, and a wave-envelope 
technique based on the method of variation of parameters was developed. 
This procedure can be used as the basis of the examination of the second 
aspect of the problem, the development of nonlinear models for the near- 
sonic region. 

The concept of sound reduction by choked inlets has been investi- 
gated experimentally at great length. Several experimental investigators 
tested actual jet engines with various shapes of centerbodies as well as 
experimental ducts to choke the flow. Surveys of the concept of the 
choked inlet were given by Lumsdaine^ and Klujber, et al^. An updated 
survey is presented here. 
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Sobell and Welliver tested a Bristol Olympus 6 jet engine; the 
inlet of this engine was choked by using a sonic block silencer. Non- 
inlet noise radiation associated with this experiment may be a reason 
why only a 12 dB noise reduction was observed. Greatrex conducted an 
experiment on an Avon engine with a bullet-shaped centerbody to choke 
the flow. He reported a 20 dB noise reduction. 

To test the effect of choking the flow on the reduction of inlet 
noise, Sawhill^ tested an ST model inlet with a translating centerbody 
and reported a 33 dB noise reduction when the throat Mach number of the 

o 

inlet was increased from 0.63 to 0.9. Cawthorn, Morris, and Hayes 
tested an SST inlet with a Viper 8 turbojet engine and a translating 
centerbody to choke the flow. They found that choking the flow resulted 
only in a 3 dB noise reduction. Using two centerbodies of different 

Q 

sizes to choke the flow of an SST inlet, Anderson obtained a 20 dB 
noise reduction at a throat Mach number of 0.77. Anderson, et al^^ 
conducted a test on an airfoil grid inlet; they used two airfoils posi- 
tioned in parallel in an inlet duct. They reported a loss of 7 % in the 
inlet recovery pressure when they attained a 27 PNdB noise reduction. 

Inlet guide vanes also have a significant effect on the noise 

reduction. Chestnutt and Stewart^ ^ conducted an experiment by using an 

accelerating inlet. They reported noise reductions up to 25 dB, due to 

the elimination of multiple pure tones, when the inlet approached 

choking conditions. The only drawback is that the noise reduction was 

accompanied by a significant reduction in the compressor efficiency. To 

determine the effect of the shape of the guide vane on the noise reduc- 
12 

tion, Chestnutt tested uncambered and tapered inlet guide vanes. He 
obtained noise reductions of about 28 dB and 36 dB for the uncambered 
and tapered guide vanes, respecti vely. Anderson, et al^^ tested radial 
vane inlets and showed a 22.5 PNdB noise reduction with a 7 % loss in 
recovery pressure. 

1 3 

Hawking and Lawson reported a large reduction in acoustic energy 
for a waisted geometry. They suggested that this reduction is due to an 
increase in the axial Mach number. Benzakin, Kazin, and Saveli 



conducted an experiment on a lined accelerating inlet. They concluded 
that the noise increases with increasing Mach number until throat Mach 
numbers of 0.6, then the noise level goes down with further increases in 
throat Mach number. 

It is clear from the above experiments that inlet choking may be an 
effective noise suppression mechanism. The amount of noise reduction 
depends on how the choking is achieved. However, the choking may be 
accompanied by a loss in the compressor efficiency. Thus, the optimum 
choking configuration is the one accompanied by no loss or a minor loss 
in the compressor efficiency. 

Many investigators studied the possibility of attaining a signi- 
ficant noise reduction with a minor loss in the compressor efficiency; 
they showed that the loss in the compressor efficiency can be minimized 
by carefully designing the centerbody. 

1 R 

Klujber reported a noise reduction when a sonic inlet is used. 

This reduction occurs when the average throat Mach number increases from 

0.5 to 1.0. He reported also that more reduction of the noise can be 

attained but with a further decrease in the inlet recovery pressure. 

1 fi 

Higgins, Smith, and Wise measured a significant noise reduction 

with a moderate loss in recovery pressure by using variable cowl inlets. 
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Koch, Ciskowski, and Garzon reported a 15 dB sound level attenuation 

with a minimum loss of aerodynamic performance when operating at an 

1 8 

average Mach number of 0.79. Miller and Abbott tested experimental ly 

an inlet with a translating centerbody to choke the flow; they reported 

19 

a 20 dB noise reduction with a pressure recovery of 98.5%. Abbott 
indicated that the most efficient method to achieve aerodynamic per- 
formance and noise reduction is to use a cylindrical centerbody at 
takeoff and a bulbshaped centerbody at approach to choke the flow. He 
reported that increasing the inlet length results in a higher recovery 
pressure for a given noise reduction. Groth tested a J-85 turbojet 
engine using a translating centerbody inlet with a radial vane. He 

measured a 40 dB reduction in a fully-choked inlet while maintaining a 

21 

recovery pressure of 92.1%. Savkar and Kazin showed that a 99% re- 
covery pressure can be attained for the same amount of noise reduction 



by proper contouring of the centerbody and careful designing of the 
22 

diffuser. Miller experimentally determined how a sonic inlet can be 
designed to have a significant noise reduction with a minimum loss of 
total pressure. 

As we see above, most, but not all, of these investigations have 
noted significant reductions of the noise level when the inlet is 
choked. The geometry of the inlet, the geometry of the centerbody and 
the operating condition seem to have an effect on the acoustic as well 
as on the aerodynamic performance. Further, most of the potential noise 
reduction is achieved by operating in the partially choked state (mean 
Mach number in the throat of 0.8-0. 9). Some investigators (e.g., 

OO C 

Chestnutt and Clark and Sobel and Welliver ) report the possibility of 
substantial "leakage" through the wall boundary layers, whereas others 
(e.g. Klujber^^) report that such leakage is minor. Although the 
experimental studies have demonstrated that the choked inlet is a viable 
technique, they have not provided insight into the physical mechanisms 
that are responsible for the noise reduction or that explain the dif- 
ferences among the several experimental results. 

Several analytical as well as numerical techniques have been 
developed for the analysis of wave propagation in uniform and nonuniform 

ducts. Surveys of these techniques were made by Nayfeh, Kaiser, and 

24 25 26 

Telionis , Nayfeh , and Vaidya and Dean . In this study, only a 

short critique is presented. 

The problem of sound propagation in a uniform duct (rectangular, 

circular, etc.), with or without mean flows, for hard as well as lined 

walled ducts, has been studied extensively. A number of parametric 

studies have been done for the case of uniform ducts, showing the 

effect of each parameter on noise attenuation. A large number of papers 

24 

are cited in the review article of Nayfeh, et al , each of which 
discusses at least one of the acoustic parameters. 

The investigation of the problem of sound propagation in nonuniform 
ducts was motivated by the experimental discoveries discussed earlier in 
this introduction. These investigations are discussed below in order of 
increasing complexity of the mean flow: no flow, one-dimensional flow, 
and two-dimensional flow. 



The problem of sound propagation in a variable-area duct with no- 

27 

mean flow was discussed for horns by Webster . He considered only the 

28 

lowest propagating mode. Stevenson extended Webster's work to inves- 
tigate the propagation of various modes. He used the method of weighted 

residuals to solve the problem of wave propagation in hard-walled horns 

29 

of arbitrary shape. Eversman, Cook, and Beckemeyer extended the 

method of weighted residuals to study multimodal propagation in a non- 

uniform lined duct. 

30 

Alfredson divided the variable area duct into a finite number of 
stepped uniform ducts. Thus, a large number of stepped uniform ducts 
are needed to provide sufficient accuracy for cases with large axial 
gradients . 

31 

Nayfeh and Telionis used the method of multiple scales to analyze 

wave propagation in ducts with slowly, but arbitrarily, varying cross 

sections and wall admittance. For the case of hard-walled ducts, the 

solution of Nayfeh and Telionis is equivalent to that of Stevenson for 

slowly-varying ducts. Nayfeh and Telionis pointed out that both of the 

solutions breakdown near cut-off; they suggested using a turning point 

analysis (see 7.3.2 of reference 32) to overcome this problem. 

33 34 35 

Isakovitch , Samuels , and Salant obtained perturbation solu- 
tions for wave propagation in ducts whose rigid walls have sinusoidal 
undulations of small amplitudes. Their perturbation expansions are not 
valid near resonance conditions; that is, whenever the wave number of 
the wall undulation is approximately equal to the sum or difference of 

oc oy 

the wave numbers of any two acoustic modes. Nayfeh ’ and Nayfeh and 
38 

Kandil used the method of multiple scales to obtain an expansion valid 
near resonance. They found that neither of the modes involved in the 
resonance can propagate in the duct without exciting the other. More- 
over, they found that coupling of an upstream and a downstream mode may 

lead to their being cutoff. 

39 40 

Quinn and Baumeister and Rice developed finite-difference 
methods to study a plane wave propagating in nonuniform ducts. We note 
that a large amount of computation will be required with these purely 
numerical techniques because a large number of grid points are needed to 



provide sufficient accuracy. The axial step must be small enough to 

resolve the smallest wave length, while the transverse step must be 

small enough to resolve the highest mode. Thus, the computational time 

increases rapidly with increasing frequency and duct length. To reduce 

the computational time for plane waves in a two-dimensional duct, 

Baumeister^^ expressed the potential function t{)(x,y,t) as ip(x,y) exp[i(kx 

- wt)], where k is a properly chosen constant, such as the wavenumber in 

a hard-walled duct. Then he solved for if^(x,y) using finite differences. 

Approximating the mean flow by a quasi-one-dimensional flow, 

42 

Powell used a multiple reflection method to study the acoustic propa- 

43 

gation through variable-area ducts. Eisenberg and Kao analyzed plane 
waves in a specially-chosen variable-area duct that yields an equation 
with constant coefficients. Davis and Johnson used a forward-integra- 
tion technique to solve the acoustic equation describing the axial 

45 

variations. Huerre and Karamcheti analyzed the propagation of the 
lowest mode by using the WKB approximation, while King and Karamcheti 
developed a second-order-accurate numerical method to solve for the 
propagation through a variable-area duct by using the method of char- 
acteristics. 

Nayfeh and co-workers studied extensively the propagation of 
various acoustic modes in ducts having slowly-varying cross sections and 
carrying general mean flows. Nayfeh, Telionis, and Lekoudis^^ dis- 
cussed the acoustic propagation in lined plane ducts with varying cross 
sections and sheared mean flow. This work was extended to annular ducts 

AQ 

by Nayfeh, Kaiser, and Telionis . The effect of a compressible, sheared 

mean flow on sound transmission through a variable-area plane duct was 

studied by Nayfeh and Kaiser . Using their method, one can determine 

the transmission and attenuation of all modes including the effect of 

transverse as well as axial gradients, but the technique is limited to 

slow variations. Moreover, the expansion needs to be carried out to 

second order in order to determine reflection and intermodal coupling of 
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the acoustic signal. Nayfeh, Kaiser, Marshall and Hurst carried out 
an experimental study of sound propagating in variable-area ducts with 
and without mean flows. Hard as well as lined duct walls were used. 

The axial variations were small. The experimental data are in reason- 
able agreement with the multiple-scales solution. 



Eversman^^ developed a theory by using the method of weighted 
residuals to determine the transmission of sound in plane nonuniform 
hard-walled ducts with mean flow. He obtained equations describing the 
axial variations of the modes. To solve these equations, one needs a 
large number of axial steps, especially as the mean Mach number ap- 
proaches unity and the frequency becomes large, leading to a rapid 
decrease in the axial wavelength. 

In summary, purely numerical techniques suffer from the requirement 

of large computation times, and they have been restricted thus far to 

cases of no-mean flow. (During the period of time that the work re- 
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ported here and in an earlier report was being carried out, finite- 
element methods for ducts with compressible mean flows were developed by 
Sigman, Majjigi and Zinn and by Abrahamson ). Analytical techniques 
have only been applied thus far to simple cases of one-dimensional mean 
flows and/or plane acoustic waves and/or slowly- varying duct geometry 
and promise to become unwieldy for more general cases. Thus, the speci- 
fic analytical and computational tools that are needed for the study of 
wave propagation in ducts involving large gradients in both axial and 
transverse directions are lacking. 

In this study an acoustic theory is developed to determine the 

sound transmission and attenuation through an infinite, hard-walled or 

lined circular duct carrying compressible, sheared mean flows and 

having a variable cross section. The theory is applicable to large as 

well as small axial variations, as long as the mean flow does not 

separate. The technique is based on solving for the envelopes of the 

quasi-parallel acoustic modes that exist in the duct instead of solving 

for the actual wave. The feasibility of this technique has been demon- 
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strated by Kaiser and Nayfeh for plane ducts with no-mean flow. 

The problem is formulated in the following section, the method of 
solution is presented in Section 3, the numerical solution is described 
in Section 4, the numerical results and discussion are presented in 
Section 5, and the conclusions are presented in Section 6. 



2. PROBLEM FORMULATION 


The transmission and attenuation of sound in hard- and soft-walled 
circular inlet ducts (Figure 1) carrying viscous or inviscid high sub- 
sonic mean flows is examined. The mean Mach number in the throat is 
near sonic; thus, the axial and radial gradients of the mean flow can be 
large. The cross section of the duct varies arbitrarily with the axial 
distance. 

The symbols used in the analysis are listed in Appendix A, except 

for a few which are used only where they are defined. All symbols are 

nondimensional unless specifically noted otherwise. Velocities, lengths, 

and time are made dimensionless by using the reference speed of sound 

c* (value in mean flow at x = 0 and r = 0), the radius R^f of the duct 
a 

in the uniform region (Figure l), and Rt/c*, respectively . The pressure 

a 

p is made dimensionless by using the density p and temperature T 

are made dimensionless by using their corresponding reference values, 
while the viscosity y and the thermal conductivity k are made dimension- 
less by using their wall values in the uniform section. In terms of 
these dimensionless variables, the equations which describe the unsteady 
viscous flow in a duct are (see for example, Schlichting ). 
conservation of mass 

If + V • (pv) = 0 (1) 

conservation of momentum 

+ V • Vv) = - Vp + V • T (2) 

conservation of energy 

p(|y V • VT) - V ’ 

= V • (kVT) + (y-I)*] (3) 

equation of state 


For a perfect gas, 
YP = pT 


( 4 ) 



For a Newtonian fluid, the dimensionless viscous-stress tensor^ and 
the dimensionless dissipation function ^ are related to v by 
= y[Vv + (Vv)^] + AV • V 

$ = T : V V = J I T. .9V ./8X . 

j=l i=l ^ ^ 

"T 

where (Vv) denotes the transpose of Vv. 

In general, the ducts carry a high subsonic, steady, sheared mean 
flow that satisfies equations (1) through (4). The presence of sound in 
the ducts results in the perturbation of the flow quantities so that 
q(r,t) = qo(r) + qi(r,t) (5) 

where q stands for any flow quantity, qo is the mean-flow part, and qi 
is the acoustic part. Substituting equation (5) into equations (1) 
through (4) and eliminating the mean-flow quantities, one obtains the 
following acoustic equations: 


1^ + V • (poVo + PiVo) = NL 

Po(gY^+ Vo *VVi + Vi • Wo) + PlVo • Wo = 

- vpi+ V • Ti + NL 


( 6 ) 

(7) 



+ Vo * VTi + Vi *VTo) + PiVo • VTo - 

Vo • Vpi + Vi • Vpo) = -J— [-^ V • (koVTi + KiVTo) 


+ (y - l)^i] + NL 

Pi - Pi I Ti 
Po Po To 


(8) 

(9) 


where and Oi are linear in the acoustic quantities and NL stands for 
the nonlinear terms in the acoustic quantities. 

No solution to equations (6) through (9) subject to general initial 
and boundary conditions is available yet. To determine solutions for 
the propagation of sound in ducts, researchers have used simplifying 
assumptions. Here, the nonlinear and viscous terms in the acoustic 



equations are neglected, and the mean flow is taken to be a function of 
the axial and radial coordinates only. Further, we neglect swirling 
mean flows. The assumption of linearization is not valid for high 
sound-pressure levels. The effects of the nonlinear acoustic properties 
of the lining material become significant when the sound-pressure level 
exceeds about 130 dB (re 0.0002 dyne/cm^), while the effects of the gas 
nonlinearity become significant when the sound pressure level exceeds 
about 160 dB. In particular, the nonlinearity of the gas must be in- 

eluded when the mean flow is transonic (i.e., near the throat ’ ). 
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Nayfeh showed that the viscous terms in the acoustic equations produce 
an effective admittance at the wall that leads to small dispersion and 
attenuation. For lined ducts, this admittance produced by the acoustic 
boundary layer may be neglected, but it cannot be neglected for hard- 
walled ducts, as demonstrated analytically and experimentally by Pestorius 
and Blackstock 

A cylindrical coordinate system (r,9,x) is introduced as shown in 
Figure 1. Since there is no swirling flow, each flow quantity qi(r,x,G,t) 
can be expressed, for sinusoidal time variations, as 


qi(r,x,0,t) = I qi (r,x)exp[-i(oot - me)] 


( 10 ) 


m=0 


where w is the dimensionless frequency. Using the above assumptions, 
one can rewrite ecuations (6) through (9) in cylindrical coordinates as 

- iwpi + -|^ (poUi + Uopi) + ^ fr ^^oPi) = 0 


Po[-iLOUi + (UoUi) + Vq |~ + Vi + pi[Up + 


( 11 ) 


+ V 


du 


° ar 
a 


L] = _ IEl 


ax 


Po[-i.v. . Ip (v.vJ . u„ . u. 1^3 . p,[v, - U„ 1^3 

ar 


av, 


av, 


Po 


r • . awi . VnWi . awi 1 im _ 

[-ICOWi + Vo ^ + Uo = - — Pl 


ar 


ax 


( 12 ) 

ax 

(13) 

(14) 
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poC-iojTi + Vo 


9x 


9r 


+ 


+ Vi 


3Tn 

9r 


+ Ui 


9To 

9x 


■] + Pi[Vo 



+ Uo 


Hjli 

9x 


(Y-DC-iwPi 


+ 


Uo 


IRl + 


9x 


Vo 


^ + 
9r 


Ui 


3Po 

9x 


+ Vi 1^] = 0 

£1. = ^ + J-L 

T f 


(15) 

(16) 


where u,v, and w are the velocities in the axial, radial, and azimuthal 
directions, respectively, and the subscript m has been suppressed. 

To complete the problem formulation, one needs to specify the 
initial and boundary conditions. The boundary conditions are based on 
the assumption that the duct wall is lined with a point- reacting acoustic 
material whose specific acoustic admittance 3 may vary along the duct. 

For no-slip mean flows, a requirement of continuity of the particle 
displacement gives 

vi - R'ui = — Pi /I + r'2 at r = R (17) 

^w^w 

where R' is the slope of the wall and the subscript w refers to values 

at the wall. The final results are desired in the form of transmission 
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and reflection matrices for a given duct section; thus the initial 
conditions consist of the successive input of each acoustic mode at the 
duct entrance. 


11 


3. METHOD OF SOLUTION 


3.1 Critique of the Existing Methods 

Since there is no exact solution available yet for equations (11) 

through (17) in ducts of varying cross sections, a number of approaches 

have been developed to determine approximate solutions to this prob- 
24 25 

lem * . These approaches include quasi-one-dimensional approxima- 

tions, solutions for slowly-varying cross sections, solutions for weak 
wall undulations, variational methods, approximation of the duct by a 
series of stepped uniform cross sections, and finite-difference and 
finite-element methods. A short critique of these approaches is in- 
cluded next (more detailed critiques are given in references 24 and 25), 
and it is followed by the proposed acoustic wave-envelope technique. 

In the quasi-one-dimensional approach, one can determine only the 
lowest mode 'in ducts with slowly- varying cross sections and cannot 
account for transverse mean-flow gradients or large wall admittance. 

In the slowly-varying cross-section approach, one can determine the 
transmission and attenuation for all modes including the effects of 
transverse as well as axial gradients, but the technique is limited to 
slow variations and the expansion needs to be carried out to second 
order to determine reflections of the acoustic signal. 

In the weak wall-undulation approach, one assumes that the dimen- 
sionless duct radius is described by R = 1 + eRi(x), where c is small. 
Thus, in this approach one can account for all effects except large 
axial variations. 

In the variational approach, one uses either the Rayleigh-Ri tz 
procedure, which requires knowledge of the Lagrangian describing the 
problem, or the Galerkin procedure (the method of weighted residuals). 
Since the Lagrangian is not known yet for the general problem, the 
Galerkin procedure is the only applicable technique at this time. 
According to this procedure, one chooses basis functions (usually the 
mode shapes of a quasi-parallel problem) and represents the pressure, 
for example, as 



CO 


(18) 


Pi = I P„(x)')’n(i".x) 
n=l 

where the are the basis functions which, in general, do not satisfy 
the boundary conditions. On expanding all flow variables in the form of 
equation (18), substituting the result into equations (11) through (17), 
and using the Galerkin procedure to minimize the error, one obtains 
differential equations describing the p^. Since the do not satisfy 
the equations and the boundary conditions, a large number of terms are 
needed to satisfy the equations and the boundary conditions and hence 
represent the solution for large cross-sectional variations; this leads 
to serious convergence questions. These problems can be minimized by 
choosing the to be the quasi-parallel mode shapes corresponding to 
the propagation coefficients k^. The functions Pp(x) vary rapidly even 
for a uniform duct: P^(x) exp(ik^x), and k^ can be very large for high 
frequency, low-order modes. Thus, small axial steps must be used in the 
computations, resulting in a large computation time, which increases 
very rapidly with axial distance and sound frequency. 

In approximating a duct with a continuously-varying cross-sectional 
area by a series of stepped uniform ducts, -a large number of uniform 
segments are needed to provide sufficient accuracy for the solution when 
the axial gradients are large. Thus, this approach is impractical in 
the present problem because an enormous amount of computation time is 
needed even for the case of a moderate number of uniform segments. 

This short discussion shows that the above techniques would cer- 
tainly either fail to produce sufficient accuracy for the present 
problem or would require large computation times. Thus, alternate tech- 
niques must be developed. The procedure to be developed has the further 
requirement that it be capable of calculating the transmission and 

reflection coefficients of the duct modes in order that the results be 

59 

compatable with the general approach developed by Zorumski . The 
generation of transmission and reflection coefficients require repeti- 
tive calculations through the duct as each mode is considered succes- 
sively as input to the duct. For each such calculation, it is necessary 
to use very small axial and radial steps to represent the rapidly- 
varying mode shapes and the axial oscillations of each mode. (In fact. 


13 



a computational difficulty exists in calculating the higher-order Bessel 
functions that represent the mode shapes in a uniform duct carrying 
uniform mean flow unless asymptotic expansions are used.) The axial 
step must be much smaller than the wavelength of the lowest mode in 
order to be able to determine the axial variations. These small steps 
could cause the error in the numerical solution to increase very rapidly 
with axial distance. Thus purely numerical methods would be impractical 
because of the excessive amount of computation time. Similar problems 
have been encountered by astronomers who developed what is usually 
called the special perturbation method in which one solves only for the 
wave envelope instead of solving for the wave itself. Here we use this 
idea to develop a wave-envelope technique for solving the present problem, 

3.2 Form of Solution 


With this approach, one uses the method of variation of parameters 
to change the dependent variables from the fast-varying variables to 
others that vary slowly. Moreover, the solution is approximated by a 
finite sum of the quasi-parallel-duct eigenfunctions. 

Thus, we seek an approximate solution to equations (11) through 
(17) in the form 


Pi 


j A^(x)4jP(r,x)exp(ijk^(x)dx) + A^(x)i|;P(r,x) x 


exp(ilk^(x)dx) 

I { A (x)i|;JJ(r,x)exp(ilk^(x)dx) + A (x)i|;JJ(r,x) x 
n=l I . / 


(19) 




exp(il k„(x)dx) 


( 20 ) 


with similar expressions for Vi, Wi , Ti, and pi, where the tilde refers 
to upstream propagation, the ^|j^(r,x) are the quasi-parallel mode shapes 
corresponding to the quasi-parallel propagation coefficients k^(x), and 
the A^(x) are complex functions whose moduli and arguments represent, in 
some sense, the amplitudes and the phases of the (m,n) modes. The 



I 

circumferential mode number m is assumed to be specified and the corres- 
ponding subscript on A, i|j, and k is not explicitly stated; each variable 
is expressed as a summation over a finite number of radial modes N, with 
n = 1 denoting the fundamental radial mode rather than the conventional 
n = 0. Since is complex, the exponential factor contains an estimate 
of the attenuation of the (m,n) mode. Thus, the envelope of the (m,n) 
mode is given by 

|A^(x)|exp[-Ja^(x)dx] 

where a is the imaginary part of k„. 

Since the i|j^ are the quasi-parallel mode shapes, they are the solu- 
tions of the following problem: 


where 
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r 
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Equation (21) - (28) can be combined to yield the following problem 
for 

1^ + [1 + ^ + 2kul3 + [f^ - - 7^]')''' = 0 (29) 

- i -rf- / = 0 at r = R (30) 

dr J ' ^ 

w 
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At each axial location, the solution of equations (29) and (30) yields 
i|;P(r;x) and its corresponding propagation coefficient kp(x). Since the 
basis functions i|^^(r;x) vary in the axial direction, they must be 
normalized in some manner to provide significance to the axial varia- 
tions of the mode amplitudes. The normalization used in this study is 

59 

the same as that defined by Zorumski : 

R 

J (r;x)]^dr = 1 
0 

Then, equations (21) - (26) are used to express the mode shapes of the 
other flow variables in terms of and k^. 


3.3 Constraints 

Since the transverse dependence in the assumed solution, equations 
(19) and (20), is chosen a priori, it cannot satisfy equations (11) - 
(17) exactly. Thus, the assumed solution must be subjected to con- 
straints. Rather than using the usual method of weighted residuals 
which forces the residuals in each of the basic equations (11) - (16) 
and the boundary condition (17) to be orthogonal to some a priori chosen 
functions, we require the deviations from the quasi-parallel solution to 
be orthogonal to every solution of the adjoint quasi-parallel problem. 

This approach assures the recovery of the results of the method of 
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multiple scales when the axial variations are slow . 

To enforce the contraints, one must define the problem adjoint to 
the quasi-parallel problem. To this end, one can multiply equations 
(21) - (26) by the functions <|)i, cj) 2 , 4 > 3 > 4>s» and (pe, respectively, 

where the (})^(r,x) are solutions of the adjoint problem, add the result- 
ing equations, integrate the result by parts from r = 0 to r = R thereby 
transferring the r-derivatives from the ip's to the cj)'s, and obtain 


R. R 

f iJj^[-itJ0(i)i - (|)6/po]dr + f ipoii^^[- 0ii(t)2 ^ k(j)i]dr + 

R 


3u, 

3r 




'r 

+ 7 (}>i]dr + c|)4 + i(y-l)aic|)5 

^ ^ X ^ V D R 

<l)6/Po]dr + / Pol|j [-1a3Cf)5 - Po(j)6]cl>" [po^ (j>l = 0 

■^0 (31) 


3r 


dr + J' 1 


I Po^'^C- 


l03(|); 


Poll) [" 00^1 
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Then, the adjoint equations are obtained by setting each of the bracketed 
terms in the integrands of equation (31) to zero; that is 


io3cf)i + 4>6 /po = 0 (32) 

/\ 

- k(j>i = 0 (33) 

-io)c{)3 + -|^ (f)2 - r 1^ (^) + 05 = 0 (34) 

~ W04 + ^ 01 =0 (35) 

ik02 - 04 + 1 (y-1)w05 + 0g/Po = 0 (36) 

Poico05 + 06 /To = 0 (37) 

Equation (31) is reduced to 

(PO0^01 + 0*^03)^^q = (pO0'^0 + 0*^03),„=[^ 


From equations (32) - (35) and (37), one can express each of the 0^ as 
a functi on of 0i ; 


02 

- ^ 

(39) 

03 

_ irTo 3 /(i)0i\ 

3r VTo^ 

(40) 

04 

_ m0i 

rw 

(41) 

05 

= ii. 
To 

(42) 

06 

= - ipoO30l 

(43) 


Using equations (39)-(43) 
i ng equation for 0i : 

X XO-i + n 

r 3r ar-* 


in equation (36), one then obtains the govern- 


Tok^ 


T om^ 


■]q = 0 


(44) 


where 

n 


±1^ 

rTo 


(45) 


It can be shown easily that equations (29) and (44) are the same; 
thus n and 0*^ satisfy the same differential equation. The boundary 
conditions on n are obtained from equation (38) by substitution for \}j^ 
from equation (23), for 03 from equation (40), and for 0i from equation 
(45). The result is 



irTo r . ,P 9rin _ irTo r . ,P 9m 
^ n + / - -^ t- ^ 1 

If one requires that q be bounded at r = 0, just as is, then the 
right-hand side vanishes. Using equation (30) to eliminate 3i|j*^/8r 
and noting that is arbitrary at r = R, one obtains 


|a-^n = 0 atr = R 
w 


(47) 


Since the boundary condition (47) is the same as the boundary condition 
(30), one can set n = without loss of generality, and hence one does 
not need to solve the adjoint problem. One needs only to solve the 
quasi-parallel problem to determine and then determine from 
rToipP 

according to equation (45). The remaining 4)'s are then determined from 
equations (39)-(43). 

Once the adjoint functions are known, the constraint conditions are 
determined as follows. On multiplying equations (11)-(16) by 
respectively, adding the resulting equations, integrating the result by 
parts from r = 0 to r = R to transfer the r-derivatives tc the <j)'s, and 
using equations (32)-(37) and (17), one obtains the following constraint: 

/" ( a a ^in 

I Hip[-iuok^Pi - ik^poUi + ^ (poUi + UoPi)] - rvoPi ~ (— ) 

+ cf) 2 ^[- iUok^poUi - ik^Pi + Po + Pl(Uo 1^ Vo|^) 


+ 1^] - Lii I;: (poVoc|)2„) + 4>3„[-iUok„PoVi + 


ar 


PoU, 


ax 


+ PoUi |^+ Pl(Vo |^+ Uo 1^)1 - VoVi 1^ (pod)3^) 

+ k^PoUoWi + 1^] _ (poVot})4„) 


'n' 


n^ 


ar 


+ <}55^[-iiJokpPoTi + (y-OiUok^pi + poUo |^ + PoUi + Pi(vo 

- u. 1^) - u. V, 1^)] - T, 


+ (y-1)Pi ^ (Vo<j)5^) 


dr + po(})i[R'ui + Pi{/r+\V^ -1)] = 0 

r:=R 

(49) 


^w^w 
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3.4 Equations Describing the Wave Envelopes 


Substitution of the assumed solution, equations (19) and (20), into 
equation (49) yields the following 2N equations for the A's; 


2N 

^mn dx 
n=l 


dA 2N 

- = J, 

n=1 


(50) 


where 


■ ‘I 


Ul,„(P0l|^U + Uoljj^) + Cj)2„(poUolpi! + 1|^E) + <|)3^(P0U0^ 




’mn 



n 

(51) 


n' "•'m' 

+ 4>4„PoUo^JJ + 4>5p^(poUoi^J - (y-l)uo^|J^)}dr]e 

- («o<) ^ . V. 1^) . gl] 

- 1^" (P0V0<t2^) + <1.3 JPOUO ^ 

+ /(vo + Uo 1^)] - Vo/ I- (po<t>3 ) + <l >0 

+ PoUo^] - ^ (PoVo<f.O^) + <(.5„[PoU„ + P„1.“ ^ 

If* 


W 


- |r (poV„<l>5^) + (y-1)*P |t (v„<t.5„) + <|.i^i(k^ - kj X 
(poi/jJJ + uoi|j^) + 4>2^i(k^ - k^)(^P + PoUoil^J^) + ^3j^^i(k^ 

- gPoUo< + K^i(k^ - gPoUo< + c^s^Kk^ - k^)uo X 

(poi/^I - (y-l)^p)| dr + po4>i^[R't|^1! + 77^ (/ 1 + R' ^ 


- 1)] 


r=R 


m^ 


n p, c ■ n 
w w 


i/k^dx 


(52) 


For convenience the upstream modes are now denoted by A^, n = N + 1, 
..., 2N, i.e.. 


^N+n ^n 


n = 1,2,3,. ..,N 
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4. Numerical Solution 


The procedures which implement the solution of Section 3 are 
described here, along with the form in which the results are presented. 

4.1 Mean-flow Model 


To evaluate the coefficients f and of the set of ordinary- 
differential equations for the mode amplitudes A^, one must specify all 
mean-flow variables, Uo ,Po >To ,Po »V q and their first partial derivatives 
with respect to both x and r. The form of g^^ as given in Eq. (52) 
would also require some second derivatives; however, those terms in- 
volving second derivatives have been integrated by parts and the re- 
sulting expression for g, which requires only first derivatives of the 
mean flow, has been used in the numerical development. Any suitable 
method of solution for the mean flow can be used in conjunction with the 
acoustic program provided it is capable of supplying the variables and 
their first derivatives and provided it supplies velocity fields which 
satisfy the no-slip boundary condition at the wall. 

For the present study, a simplified model of the mean flow has been 
employed. This model uses one-dimensional gas-dynamics theory to des- 
cribe the mean-flow variables in the inviscid core; thus Uo,PosTo and po 
are constant across the duct cross section except in the region of the 
wall boundary layer. The program permits one to select one of two 
options: the radial velocity Vo can be set equal to zero, consistent 
with the one-dimensional theory, or the radial velocity can be calcu- 
lated to be a linear function of r, consistent with the mean-continuity 
equation and the flow tangency condition at the wall. The velocity 
profile in the wall boundary layer is taken to be a quarter-sine pro- 
file; that is. 


= sin[7r(R - r)/26] r _> R - 5 
^ = 1 r < R - 6 


The temperature profile is related to the velocity profile by 

T - T . 

f] _ iLo.1 

- - _ ^ u 

c c c c 


^ = 1 + ^M^[l - (^)2] + 

c ‘ 

ad/Tc = 1 ^ 


(53) 
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where the subscript c refers to values in the inviscid core, is the 
wall temperature, is the adiabatic wall temperature, 6 is the bound- 
ary-layer thickness, ri is the recovery factor, and y = is the ratio 
of the gas specific heats. The quarter-sine profile is a good approxi- 
mation for acoustic studies for cases of laminar mean flow^^. 

The temperature-profile formula, Eq. (54), is a good approximation 

for laminar boundary layers in zero pressure gradients; in this case, 

1 / 

the recovery factor ri is taken to be equal to Pr for best results, 

although many boundary- layer studies have been based on ri = 1 . Since 

the present study is concerned with variable-area ducts, the use of Eq. 

(54) should be considered a rough approximation only, providing the user 

with a means of checking the sensitivity of the results to changes in 

the mean- temperature profile. To this end, three options have been 

provided in the program for selection of the wall temperature: the wall 

temperature can be set to any specified constant value, it can be set 

equal to the inviscid core value, T^, or it can be set equal to the 

adiabatic wall temperature, ^ad- It is noted that if one selects the 

second option, T = T , and also sets ri = 0 (no aerodynamic heating) 
w c 

the mean temperature will be constant across the duct width and the 
temperature-profile refractive effect will be eliminated; this is of 
value for attempting comparisons with one-dimensional results or in- 
compressible mean-flow results. 

The axial variation of the boundary-layer displacement thickness 6 i 
is assumed to be known and is specified in the program by a simple 
polynomial variation: 


61/610 = [1 + bi(x/L) + b 2 (x/L) 2 ]R (55) 

The displacement thickness and the Mach number within the uniform core 
have prescribed values, 610 and M^ , at x = 0 ; the subsequent axial 

C 0 

variations of 6 and are calculated within the program from the defi- 
nition of displacement thickness and from mass-flow considerations. The 
one-dimensional gas-dynamics theory provides the axial variation of T^, 
p^, u^, etc. and the boundary- layer profiles are computed from Eqs. (53) 
and (54). Thus, the mean flow within the duct is prescribed completely 

by input of the values of M , 61 0 >bi ,b 2 ,ri , y» and T, , to the program. 

Co w 



4.2 Wall Admittance 


The wall of the duct is assumed to be lined and two options are 
available to the user: the duct has a point-reacting liner of constant 

properties whose specific admittance is described by 

e = '' — ^ — (56) 

R_(l-i(jo/o3o ) + i cot(ojd/T / ^) 

G W 

where is the resistance of a facing sheet, ojq is the characteristic 
frequency of the facing sheet, and d is the depth of the backing cavi- 
ties in the liner; or the admittance of the liner can be taken to vary 
from a specified value Bo at x = 0 to a specified value B|_ at x = L 
according to the expression 

3 = Bo + (3 l - So)(3 - 2x/L)(x/L)" (57) 

This latter expression gives a continuously varying admittance with 
dB/dx = 0 at both x = 0 and x = L. 

4.3 Parallel-Duct Eigenfunctions 

To calculate the changes in the amplitude of the acoustic wave 
first requires the eigenfunctions The number of radial modes to be 
considered and the mode eigenvalues at x = 0 must be specified as input 
to the program. The quasi-parallel acoustic equations (29) and (30) are 
solved at each axial position by using a fourth-order Runge-Kutta 
forward-integration technique and by employing a Newton-Raphson pro- 
cedure to determine the eigenvalues k^. The accuracy of the Runge-Kutta 
procedure should assure good numerical results with a minimum of grid 
points across the duct width. 

To determine the coefficients g^^ of equation (52), one has to 
evaluate the axial gradients of the wavenumber, k, and of the eigen- 
functions These axial derivatives can be obtained by differentiat- 
ing equations (29) and (30) with respect to the axial co-ordinate x, or 
by using a simple finite-difference quotient such as 

dk^ ~ *^x+Ax ~ *^x-Ax 

2A3; (58) 



Although the first approach is more elegant and inherently more accurate 
it is somewhat tedious to implement for this case and would lead to 
longer computation times. Thus the second approach is used even though 
it introduces truncation errors of order Ax^; convergence checks using 
calculations with different values of Ax have shown that this is not a 
problem. 

4.4 Numerical Integration 

The adjointpfunctions are found from the quasi-parallel-flow 
variables -||— , and k by using the relations (39)-(43). The coef- 
ficients f|^^ and g^^^ are then evaluated from equations (51) and (52). 

The integrals across the duct in these expressions are evaluated using 
Simpson's rule which yields a high level of accuracy. The axial inte- 
grals /k^dx are evaluated with the trapezoid rule which is of sufficient 
accuracy to obtain good convergence with changes in Ax. 

Writing equation (50) in matrix form, FdA/dx = GA, and solving for 
dA/dx, one obtains 

^ = F"''gA (59) 

where A is a column vector whose elements are the A^. A Runge-Kutta 
forward- integration technique is used to solve equations (59) for the 
function A at each axial station. Since the problem is linear, one can 
determine the solution for any problem subject to general boundary 
conditions at the two ends of the duct by a linear combination of 2N 
linearly independent solutions. 

The linearly independent solutions are obtained by setting all mode 
amplitudes except one equal to zero at x = 0 and integrating equation 
(59) to x = L, One such integration for each of the 2N modes allows one 
to obtain the transfer matrices TRi, TRa, TR3, TRi+ which are defined by 
B'^(L) = TRiB'*'( 0) + TR2B"(0) (60a) 

B"(L) = TR3B'*‘(0) + TR4B“(0) (60b) 

where B (x) is a column vector of the amplitudes A^e of the 

right-running modes and B"(x) is a column vector of the amplitudes 
~ i/k dx 

A e of the left-running modes. The transfer matrices thus allow 



one to calculate the complex mode amplitudes at x = L from those at x = 
0. Following reference 59, we derive results in the form of trans- 
mission and reflection coefficients for the variable-area segment being 
considered. The transmission and reflection coefficients relate the 
complex magnitudes of the outgoing modes to those of the incoming modes, 

B'^(L) = t'-’V{ 0) + R*-’4‘(L) (61a) 

B'(0) = T°’‘-B'(L) + R°’V(0) (61b) 
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and are calculated from the transfer matrices by 


rO,L 

^0,0 

)L,L 


-L,0 


= - TR'^TR3 
= TR^TR’^ 

= TR^ + TR^R^’^ 


(62) 


The reflection coefficients are the negative of those defined in 
reference 59 as a consequence of the use of the positive sign on the 
term in equation (19). The (m,n) term of T^’ represents the trans- 
mission of the m ^ radial mode at x = L due to the n radial mode being 
incident at x = 0, etc. The requirement that the procedure be able to 
calculate these transmission and reflection coefficients makes a direct 
numerical procedure undesirable for this study. The wave-envelope, 
eigenfunction-expansion procedure developed in Section 3 is better 
suited to the the necessary repetitive calculations with each mode as 
input than a direct numerical approach would be. Further, it is noted 
that the transmission and reflection coefficients are general; no 
assumption about the nature of the source input to the duct has been 
made. 


4.5 Acoustic Pressure Profiles 

Although the procedure developed here is primarily intended to 
obtain transmission and reflection coefficients, the acoustic pressure 
distributions can be constructed if the input to the duct section is 
specified. Specifically, one must designate the values of B (0) and B 
(L) which are the amplitudes of the right- running modes that are in- 


cident on the duct section at x = 0 and the amplitudes of the left- 
running modes that are incident on the duct section at x L, respecti- 
vely. Equation (61b) then yields the value of B“(0). A generalization 
of equation (60) , 

B'^(x) = TR^B'^(O) + TR2B"(0) 

B’(x) = TR3B‘*‘(0) + TR^B"(0) 


(63) 


where the transfer matrices are functions of the axial coordinate x, is 
used to calculate the mode amplitudes throughout the duct. The program 
has been constructed to store the necessary values of the transfer 
matrices. Finally, the pressure distribution across the duct at each 
axial position is obtained from equations (10) and (19) and the defini- 
tion of B^ and B”: 


Pi(r,x,6,t) 


N ^ ^ N 

I B^(x)i()P(r,x) + I 
n=l " " n=l 


B-(x)iP(r,x)le’('"®-“t) 

(64) 


where the bracketed terms describe the spatial distribution of interest. 
A similar expression is evaluated for the acoustic particle velocity 
(axial component) and the moduli and arguments of these quantities are 
output of the program. 

The accuracy of these acoustic pressure and velocity profiles de- 
pends on the accuracy of the eigenfunction-expansion procedure that 

forms the basis of this study. General guidelines are available from 
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parallel-duct studies, for example, Hersh and Catton and Unruh and 

CO 

Eversman ; eigenvalues and attenuation rates (transmission and reflec- 
tion coefficients in the variable-area case) converge more rapidly with 
an increasing number of basis functions than do the eigenfunctions 
(pressure profiles in the variable-area case). By using the parallel- 
duct eigenfunctions as the basis functions, we should minimize any 
convergence problems if all cut-on modes are included in the analysis; 
the magnitudes of any cut-off modes should be sufficiently small so as 
not to cause a large error. In addition, the integrabil ity constraint 
developed in Section 3.3 assures that the governing equations are 
satisfied "on the average" at each cross section; that is, the average 



error at any cross section will be zero even if an inadequate number of 
modes is used. Thus the overall characteristics of the duct such as 
transmission and reflection coefficients are more accurate than the 
details of the acoustic pressure distributions. 

4.6. Energy Flux 

When the acoustic pressure and velocity distributions are obtained, 
an evaluation of an energy-flux expression is possible. We have used 
the expression developed by Morfey^^ (and others) in the computer pro- 
gram developed for this study; the reader is cautioned, however, that in 
the presence of a mean flow with strong vorticity, the "acoustic energy" 
as defined by this expression is not conserved. Morfey points out that 
there are production terms present unless the mean flow is irrotational . 

Alternative expressions for the energy flux are also of questionable 
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value for the current study. The expression developed by Eversman 

(and others) is restricted to uniform mean flows and the expression 

fifi 

developed by Mbhring while permitting rotational mean flows does not 
allow for mean flows in which vorticity is developed as the mean flow 
moves in the axial direction and thus is unsuitable for the variable- 
area cases of interest in this study. Whether any of these expressions 
are of value for studies of variable-area ducts with developing mean 
flows that have wall shear layers is questionable. For a few simple 
cases with no mean flow or low-speed mean flow any of the expressions 
provides a check on the validity of the method of solution and its 
computational implementation. 
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5. Sample Cases 


The computer program described in the previous section has been 
developed and used to provide a number of sample results. These sample 
results serve three purposes: to test the validity of the program and 

the method of analysis; to indicate the range of mean-flow and acoustic 
parameters for which the analysis is useful and simultaneously mark 
potential problems that may develop; and to provide a set of test cases 
that can be used by anyone who wishes to become familiar with the use of 
the program. To assist with this last requirement, each case contains a 
complete set of the input variables that were used to run the case. 

All cases were run in single precision on the IBM 370 at VPI&SU; 
many of them have also been checked with a version of the program 
implemented on the CDC Cyber 175 at Langley Research Center. In most of 
the test cases, the duct radius is assumed to have a simple converging- 
diverging variation with axial distance: 

R = 1 + azC-l + cos(2Trx/L)] (65) 

where a 2 specifies the magnitude of the variation in the outer wall and 
L is the length of the duct. If one wishes to consider only a con- 
verging (or diverging) duct, the program allows the calculations to be 
terminated at x = L/2. For other geometric variations, one must alter 
the geometry subroutine in the program; these alterations are simply 
accomplished, and one set of cases included here requires such changes. 

5.1 Definition of Input Variables for the Program 

To permit the reproduction of the results contained in subsequent 
sections, the Fortran variables used as program input and, when appro- 
priate, their correspondence to the variables used in the analysis are 
defined here: 


Variables that de fine the mean flow 

Fortran Variable 
XMCO 
GAM 


Definition 

M^ , Mach number at the duct entrance 
Co 

Y, Ratio of specific heats 
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RCF 

DISP20 

B2,B3 

CWT 

NTEMP 


NRADVL 


A2 

XF 

NCUT 


RE 

WO 

DIS 

BETAO 

BETAL 


W 

M 

NM 


Pi , Recovery factor 

610 , Displacement thickness at x = 0 

bi,b 2 , Coefficients in eq. (55) describing the 

displacement thickness variation 

Constant value of the wall temperature 

= 1 , = CWT 

w 

" " ^ad 

= 3. T„ = 

= 0 for zero radial velocity 

= 1 for linear variation of the mean radial 
velocity 

Variables that define the duct properties 

32 , Coefficient that determines the magnitude 
of the variation in the duct radius 
L, the duct length 

= 1, Calculations for the full converging, 
diverging section are carried out 
= 2, Calculations for only the converging 
portion of the duct are considered 
Rg, Liner resistance 

0 ) 0 , Characteristic frequency of the liner 
d. Cavity depth of the liner 
3o» Liner admittance at x = 0; used only if 
RE = 0 

3|_, Liner admittance at x = L; used only 
if RE = 0 

Variables that define the acoustic signal 

0 ), Circular frequency 
m, Circumferential mode number 
2N, The total number of parallel-duct modes to 
be used in the calculations 



IVK The initial approximations for the eigen- 

values of the parallel-duct modes at 
X = 0; an array. 

ERR The error bound for convergence of the itera- 

tions for the parallel-duct eigenvalues 

Variables that define the numerical procedure 

NSC The number of steps across the duct in the 

uniform mean flow core; must be an even 
number 

NSB The number of steps across the duct in the 

mean- flow boundary layer; must be an 
even number 

NOS The number of axial segments into which the 

duct length is divided. Flow variables 
and acoustic profiles are printed at the 
end of each segment. 

NSS The number of axial steps that each segment is 

sub-divided into; must be either 2 or 4. 

Ax = 1/(N0S NSS) 

Variables that control output 

A number of variables may be printed at the end of each duct seg- 
ment to serve as diagnostic tools should a problem with the program 
develop. Each of the following parameters is set to zero if the corres- 
ponding printout is not desired and set to unity if it is desired. 

NRITR Prints the initial iterations on the 

parallel-duct eigenvalues at x = 0 

NREGN Prints the converged values of the parallel- 

duct eigenvalues at the end of each 
segment. 

NRAMP Prints the amplitudes of each duct mode 

at the end of each duct segment. 


Variables for calculating the acoustic pressure profiles 


NINPCS Specifies the number of cases for which 

the acoustic pressure and velocity 
distributions are to be calculated 
BPO B^(0); Mode amplitudes of the right- 

running modes at x = 0; an array 
BML B"(L); Mode amplitudes of the left- 

running modes at x = L; an array 


5.2 Symmetry Checks 


Testing of the method of solution is somewhat difficult since there 
are no available solutions for acoustic propagation through variable- 
area ducts carrying a high-speed, compressible mean flow, except for 
some one-dimensional solutions. However, a few tests can be conducted 
for simple flows to determine internal consistency of the program. 

For straight, uniform ducts, the wave-envelope amplitudes are 
correctly calculated to be constant throughout the duct, and the routine 
for solving the parallel-duct eigenfunctions and eigenvalues produces 
the expected results. In addition, a number of checks on the trans- 
mission and reflection coefficients have been performed for both straight 
and variable-area ducts with and without liners and with and without 
mean flows. A few of these symmetry checks are presented here. 


(a) Straight, lined duct with no mean flow. 


Since the duct is uniform, the reflection coefficients should be 
identically zero, and the off-diagonal terms of the transmission ma- 
trices should be zero; further, should be equal to because the 

left-running and right-running modes are identical. Results for the 
transmission and reflection coefficients for one such case are 


f ,0 ^ ^0,L ^ 
,0,0 ^ ,L,L ^ 


■-.241041 - .7819581 
-.000022 - .0000181 

■ .000000 - .OOOOOOi 
-.000003 - . 000001 i 


.000014 + .0000121' 
.551753 + .6712031. 

.000002 + . 000001 i 
.000000 + .0000001 



The expected results are achieved to a suitable degree of accuracy with 
errors in the intermodal coupling being larger than those in reflection 
of the modes. The input data for this case are W=5, M=0, NM=4, 
IVK = 5.35 + O.li, 3.58 + .07i, - 5.35 - O.li, -3.58 - .07i, ERR = 

10"^, XMCO = 0.0, GAM = 1.4, RCF = 1.0, CWT = 1.0, NTEMP = 1. DISP20 
.01, B2 = B3 = 0.0, NRADVL = 0, A2 = 0.0, XF = 2.0, NCUT = 1 , RE = 0.8, 
WO = 15, DIS = .05, NSC = 16, NSB = 2, NOS = 20, NSS = 2, NINPCS = 0. 
Note that one must specify a non-zero boundary- layer thickness which, 
for cases of no mean flow, is strictly an expedient for setting up the 
numerical steps across the duct. In this case the number of steps in 
the "boundary layer" is set to the minimum possible value since the 
"boundary-layer thickness" is small. 

(b) Variable-area lined duct with no mean flow. 


In this case, reflection and intermodal coupling in the transmis- 
sion coefficients is expected to be non-zero; however, the symmetry of 
the situation still requires T^’*" = and R^*^ = R*"’^. For one such 
case, the results are 


tL,o ^ 


-.17802 - .63795i -.33815 - .16619i 

-.50851 - .24930i .57103 - .24494i 


T 


0,L = 


-.17803 


. 63804i 


-.33818 - .16623i 


- .50856 - .24936i 


.57110 - .24497i 


Errors are typically in the fifth decimal place with the maximum error 
in the fourth decimal place. The reflection coefficients agree to five 
decimal places, and 

Q Q r. 00426 + .01666i - .03609 + .04866i 
L. 05460 + .07300i - .14420 - .14625i 


The input data for this case are the same as those in part (a) except 
that A2 = 0.1. 


(c) Variable-area lined duct with mean flow 


In this case the transmission coefficients for a case when 
M > 0 (mean flow from left to right) should be the same as the coeffi- 
cients for a case when M < 0 (mean flow from right to left) pro- 

vided all duct properties and the mean flow are symmetric about x = L/2. 
Similar arguments can be made for the reflection coefficients. An 
example of such results yields 


for M = - 0.3, 
Co 


-L,0 _ 




and for = - 0.3 
Co 

tO,L _ 


[: 


51771 + .29426i 
38645 - .31743i 


51790 + .29418i 
38650 - .31753i 
0,L 


.37739 - .27368i 
-.06760 + .36536i 


.37744 - .27382i 
-.06748 + .36548i 
-L,0 


The agreement between T"'*'" for Mo = - 0.3 and T""’"" for Mo = + 0.3, is 
comparable to that shown above and the agreement for reflection coeffi- 
cients is better. The results are relatively insensitive to changes in 
the error bound on the eigenvalues of the parallel duct; changes in the 
numerical step sizes will alter the values of the transmission coeffi- 
cients but do not change the general conclusions concerning the self 
consistency of the results. Since the largest errors are in the fourth 

decimal place, the results are considered to be satisfactory. 

-L,0 


The input data for the above results for T * are W = 10, M = 0, 

= 4, IVK = 7.2+.081, 5.3+.22i, -13.9-.04i, -12.3-.20i, ERR = 10"^, 

XMCO = 0.3, GAM = 1.4, RCF = 1.0, CWT = 1.0, NTEMP = 1, DISP20 = .01, B2 


= B3 = 0.0, NRADVL = 0, A2 = 0.12, XF = 2.0, NCUT = 1 , RE = 0.8, WO = 

15, DIS = .05, NSC = 16, NSB = 8, NOS = 20, NSS = 2, NINPCS = 0. For 
calculating T^’^ the input data was the same except for XMCO. = -0.3, 
and IVK = 13.9 + .04i, 12.3 + .20i, -7.2-.08i, -5.3-.22i. The mean flow 
Mach number in these cases reaches a value of .61 at the throat of the 
duct, sufficiently large that compressibility effects in the mean flow 
are important. 



5.3 Comparison with one-dimensional theory 

For a one-dimensional mean flow and acoustic disturbance, results 
are available in the literature; for example, Myers and Callegari^ 
use such calculations to indicate the singular behavior of linear acoustic 
theory as the mean Mach number approaches unity. For their calcula- 
tions, the duct area is described by a polynomial 

= + 8 (^- 1 }^ - 16 (^- 1 )"] 

and results are presented in the form of the ratio of the acoustic pres- 
sure amplitude to the incident pressure amplitude as a function of 
distance through the duct. 

The wave-envelope method clearly is not a one-dimensional model; 
however, as many two-dimensional effects as possible have been sup- 
pressed to effect the closest possible agreement between conditions in 
the one-dimensional calculations and the wave-envelope calculations. 

Only one mode propagating in each direction is included, and the duct 
wall is taken to be rigid. The boundary- layer displacement thickness is 
set to a very small value (.001) in an attempt to minimize any refrac- 
tive effects on the acoustic signal. In order to obtain the same axial 
variation of the mean flow in both calculations, the radius of the duct 
for the wave-envelope calculations has been taken to be 

R = ||[1 + 8(f - Iv - 16(^ - + S, 

and the displacement thickness is taken to be constant 

<5 1 - 6 1 0 “ • 001 

These expressions require that the geometry subroutine of the program be 
altered if one wishes to reproduce the results of this section. To 
suppress two-dimensional effects, the mean radial velocity has been set 
to zero; the wall temperature is taken to be equal to the value at the 
centerline and the recovery factor has been set to zero in order to 
eliminate any refractive effect from a transverse temperature gradient. 
Despite .these restrictions, the wave-envelope method is still basically 
two dimensional: for example, the acoustic particle velocity normal to 
the wall is required to be zero whereas a one-dimensional model is 
equivalent to a zero acoustic velocity normal to the duct centerline. 



Hence it is of interest to examine the differences between the results 
of the two approaches. In general, it has been found that the magnitude 
of the acoustic pressure is larger near the throat of the duct than is 
predicted by one-dimensional theory, especially when the throat Mach 
number is large. In addition, as the throat Mach number reaches high 
subsonic values, a substantial refractive effect develops, even from the 
very small boundary layer that was used in this study. An example of 
these effects is shown in Fig. 2. The converging duct section results 
in a higher acoustic pressure in the throat region than is predicted by 
one-dimensional theory. The incident signal at x = 0 propagates down- 
stream, and an upstream wave is present as a consequence of reflection 
within the variable-area duct. In the throat region the upstream wave 
has a very short wavelength and thus the boundary layer produces a 
significant refraction of the upstream wave toward the duct centerline. 
The acoustic propagation in the throat region is not a one-dimensional 
process . 

The effect of removing some of the "one-dimensional" restrictions 
has been examined. The results of Fig. 3 have been obtained with a 
constant wall temperature throughout the duct, T, = 1.0; thus the wall 
is hotter than the mean flow in the vicinity of the throat and a temper- 
ature-prof ile refractive effect occurs. This results in a further in- 
tensification of the acoustic pressure in the throat region. Finally 
the influence of a mean radial velocity has been included, and the 
results are shown in Fig. 4. The mean radial velocity component has a 
strong effect on the acoustic pressure, reducing the pressure amplitude 
at the throat and bringing the two-dimensional results into closer 
agreement with one-dimensional theory than occurred in the two previous 
cases. Input data for reproducing the case in Fig. 2 are W = 1.009, M = 
0, NM = 2, IVK = 7.77 + O.Oi, -14.38 + O.Oi, ERR = 10"^, XMCO = .299, 

GAM = 1.4, RCF = 0.0, NTEMP = 3, DISP20 = .001, NRADVL = 0, XF = 1.0, 
NCUT = 1, RE = 0.0, BETAO = 0.0 + O.Oi, BETAL = 0.0 + O.Oi, NSC = 14, 

NSB = 10, NOS = 40, NSS = 2, NINPCS = 1 , BPO = 1//2, 0.0, etc., BML = 
0.0. For Fig. 3, the data are the same except that NTEMP = 1 and CWT = 
1.0. Finally the results in Fig, 4 have the same input data as those 
for Fig. 3 except that NRADVL = 1. 



5.4 Comparison with Finite-Element Results 

A number of test cases have been set up for comparing the results 
of the wave-envelope procedure with the results from the finite-element 
procedure developed by Abrahamson . The finite-element results have 
been supplied by H. Lester of Langley Research Center. Currently, 
comparisons have been made for two of the cases, and additional ones are 
being carried out. 

The finite-element method was developed for a compressible mean 
flow, but is currently implemented with an incompressible model; hence 
the cases chosen for comparison are for low-speed mean flows or for no 
mean flow. The mean flow varies from M = -0.1 at x = 0 to M = -.158 
at the throat, and the mean boundary- layer thickness is assumed to be 
10% of the local duct radius; the wave-envelope solution is based on a 
uniform mean temperature distribution across the duct. The wave-envelope 
method specifies the mode amplitudes incident on the duct section at x 
= 0, specifies no input at x = L and solves for the acoustic pressure 
distribution throughout the duct. The pressure distribution that re- 
sults at X = 0 is used as input for the finite-element program, which 
specifies a poCo impedance at x = L to approximate the no-input boundary 
condition at that station. The pressure profiles from the two methods 
at x/L = 1/4, 1/2, 3/4 and 1, as well as the axial variation of the 
pressure at the wall, are then compared. 

The two cases currently available correspond to an input signal of 
500 Hz in a 2-meter-diameter rigid duct with the mean speed of sound 
being 344.4m/sec. These conditions yield a dimensionless circular fre- 
quency of 03 = 9.12 at which there are three cut-on modes propagating in 
each direction. The comparison for no mean flow is shown in Fig. 5. 

The incident signal at x = 0 is the fundamental, plane mode and reflec- 
tion in the duct results in the acoustic distribution shown in the left 
portion of the figure. Considerable variation of the acoustic profiles 
takes place as the signal propagates through the duct. The trends from 
the two methods of solution are the same, and the results agree reasonably 
well except at the center of the duct near the throat. The variation of 
pressure along the outer wall is also quite pronounced, and both methods 



yield the same trends: the peak magnitude occurs ahead of the minimum 
area and a local minimum in pressure occurs close to the minimum area. 
Although the agreement is not exact, overall it is very satisfactory. 

The transmission coefficient and reflection coefficient 
for this case are 


and 

r0,0 ^ 


.823-.519i 

.118-.045i 

.135-.181i 

.003-. 0121 
.011+.020i 
-.077+. 0211 


.106-. 0401 
-.703-.5321 
.121+.5141 

.010+.0181 

-.027-.0151 

.082-.094i 


.085-. 1141 
.085+. 3611 
-.342+. 7271 

-.049+. 0141 
.058-. 0671 
.276+. 1801 


The intermodal coupling coefficients (the off-diagonal terms of T^’^) 
are an order of magnitude larger than the reflection coefficients, and 
it is this coupling between modes that is primarily responsible for the 
variations of the acoustic profiles throughout the duct. Despite the 
small values of the reflection coefficients, they have a noticable 
effect on the acoustic pressure distribution as seen in Fig. 5 by 
comparing the incident signal with the total signal at x = 0. The input 
data to reproduce this result are W = 9.12, M = 0, NM = 6, IVK = 9.12, 
8.28, 5.83, -9.12, -8.28, -5.83, ERR = 10’^, XMCO = 0.0, GAM = 1.4, RCF 
= 0.0, NTEMP = 3, DISP20 = .0363, B2 = B3 = 0.0, NRADVL = 0, A2 = 0.1, 

XF = 2.0, NCUT = 1, RE = 0.0, BETAO = BETAL = 0.0, NSC = 20, NSB = 2, 

NOS = 16, NSS = 2, NINPCS = 1 , BPO = 1//?, 0.0, etc. 

A similar comparison has been made when the input signal propagates 
upstream against the mean flow. The results are shown in Fig. 6 and are 
very similar to those for no mean flow. The mean boundary layer causes 
a refraction of the disturbance away from the duct wall as expected; the 
differences between the finite-element and wave-envelope results are 
very large in the throat region near the duct centerline but otherwise 
are comparable to those found in the case without flow. Both methods 
"conserve energy" to a suitable level of accuracy. The transmission 
matrix T^’^ is 


T 


L,0 


1 -.314+. 9091 
-.050+. 1241 
-.005+. 2861 


-.045+. 1131 
.816+. 0871 
-.449-. 4571 


-.003+. 1901 
-.329-. 3351 
-.050-. 7631 

_ 


Again the intermodal coupling in transmission is quite strong. Input 
data for reproduction of the results are the same as those for the no- 
mean-flow case, except for XMCO = -0.1, NRADVL = 1, NSB = 8, and IVK = 
10.13, 9.29, 6.85, -8.29, -7.45, -5.01. 

At a frequency of 1000 Hz or oj = 18.24, there are six cut-on modes 
at the entrance, but the sixth mode cuts off within the duct. This 
situation presents difficulties for the wave-envelope calculations since 
it breaks down at a hard-wall cut-off point; it is also noted that if 
the sixth mode is not well cut on the poCo impedance used as a boundary 
condition at x = L in the finite-element program may not be adequate. 

The breakdown of the wave-envelope method at hard-wall cut off is the 
major restriction on the use of the wave-envelope procedure. The re- 
sults can be obtained with only five modes, but one would expect this to 
introduce errors. Thus, the wave-envelope procedure is better suited 
for lined ducts than it is for rigid ducts. 


5.5 Straight Ducts with Variable Liners 

Several cases have been considered for wave propagation through a 
straight duct in which the wall admittance varies from zero at x = 0 to 
= 0.6 + 0.6i at X = L. Other conditions are basically the same as 
those discussed in Section 5.4, and ultimately a comparison of these 
results with the finite-element results will be made. Two results are 
given here, for oj = 9.12 and co = 18.24, both for an incident signal at 
X = 0 that consists of two modes propagating upstream. In these cases, 
the reflection coefficients are so small that they do not affect the 
acoustic profiles. The results for the lower frequency are shown in 
Fig. 7 and those for the higher frequency in Fig. 8. The refraction of 
the signal away from the wall due to the upstream propagation is more 
pronounced in the high-frequency case. The acoustic amplitude along the 
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wall is reduced significantly as the wave propagates through the duct, 
but the reduction in energy flux is small in both cases (1.46dB and 
.34dB), The main' effect of the variable liner properties is to couple 
the modes in transmission as shown by the transmission matrix for co = 
9.12: 


tL,o ^ 


.448+. 782 i 
-.185+.132i 


.333+.052i 

.307-.632i 


-.001+.020i 


-.008-.058i 


.007+.024i 

~.158+.062i 

.434-.100i_ 


The first two modes interact strongly, whereas only small amounts of the 
third mode are generated by the presence of the first two. Input data 
for the low frequency case are W = 9.12, M = 0, NM = 6, IVK = 10.13, 
9.29, 6.85, -8.29, -7.45, -5.01, ERR = 10"^, XMCO = -0.1, GAM = 1.4, RCF 
= 0.0, NTEMP = 3, DISP20 = .0363,82 = 83 = 0.0, NRADVL = 1 , A2 = 0.0, XF 
= 2.0, NCUT = 1, RE = 0.0, BETAO = 0.0, 8ETAL = 0.6 + 0.6i, NSC = 20, 

NSB = 8, NOS = 16, NSS = 2, NINPCS = 1, 8P0 = 1/v^, 0.1-0.2i, 0, etc., 
8ML = 0.0. The higher frequency corresponds to W = 18.24, IVK = 20.27, 
19.86, 18.86, 17.17, 14.5, 9.93, -16.58, -16.17, -15.18, -13.48, -10.81, 
-6.25, with all other data the same as in the low frequency case. 


5.6 General, Axisymmetric Case 

A case is examined in which the incident signal propagates down- 
stream in a converging-diverging duct in which the mean flow reaches a 
Mach number of .465 at the throat and the boundary-layer thickness grows 
rapidly downstream of the throat. In addition, the duct is lined with 
a variable-admittance liner, changing from 3o = .5-.2i to 3|^ = .2-.li. 

It is assumed that the mean flow has reached adiabatic conditions with 

the duct walls: T = T^ .. The basic physical properties are defined by 

w ao 

a 1000 Hz signal propagating in a duct whose diameter varies over a 
60.96 cm (2‘) axial distance from 30.48 cm (!') at the entrance to 25.4 cm 
(10") at the throat. Using Co = 335.28 m/sec. (1100 ft/sec), one 
obtains oo = 2.86, a frequency at which only one cut-on mode propagates 
in each direction. The transmission and reflection coefficients for 
this case are 
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I 


= .0683-. 16361 
= -.0007+. 00021 
= -.0023-. 00221 
= -.0004+. 00011 

Substantial reduction of the acoustic signal occurs In this case as 
shown by the small magnitudes of the transmission coefficients. An 
upstream-propagating signal would be almost eliminated. For the first 
right-running mode incident at x = 0, the pressure profile at x = L also 
is quite small compared to that at x = 0, as shown in Fig. 9; the energy 
flux Is reduced by approximately 15.3dB. Input data for this case are 
W = 2.86, M = 0, NM = 2, IVK = 2.7+.481, -3.16-.76i, ERR = 10‘^, XMCO = 

0.3, NRADVL = 1, GAM = 1.4, RCF = 0.8, NTEMP = 2. DISP20 = .02, B2 
= -2. B3 = 4, A2 = .08333, XF = 8.0, NCUT = 1 , RE = 0.0, BETAO = 0.5-0.21, 
BETAL = 0.2-0.11, NSC = 16, NSB = 8, NOS = 25, NSS = 4, NINPCS = 1, BPO 
= 1//2, 0.0, etc. , BML = 0.0. 

5.7 Spinning-Mode Case 

All the previous sample calculations considered axi symmetric 
acoustic waves. The program also can be used to study cases involving 
spinning modes as illustrated by a case for m = 2. The input signal of 
frequency w = 7 propagates downstream in a variable-area duct with 
constant liner properties. The acoustic pressure profiles at several 
axial positions along with the wall pressure variation are shown in Fig. 

10; these profiles correspond to the lowest radial mode being incident 
at X = 0 and are calculated with the two cut-on radial modes propagating in 
each direction. Transmission and reflection coefficients in this calcu- 
lation are 


^L,0 ^ r. 475+. 0971 

-.030-.066i 

-.078-. 1441 

.013+.0621 

j^0,0 ^ r.008-.026i 

.108+.011i 

|_. 257+. 0311 

.336+. 5671 


It is seen that the second downstream mode is strongly reflected into 
two upstream-propagating modes with relatively little of the mode 
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passing through the duct section. This effect is not seen in Fig. 10 
since it corresponds to the case when only the lowest mode is incident 
on the duct section. Input data for these calculations are W = 7.0, M = 
2, NM = 4, IVK = 6.4+.32i, 2.47+.25i, -8.4-.63i, -4.04-.28i, ERR = 

10"^, XMCO = 0.1, GAM = 1.4, RCF = 0.0, NTEMP = 3, DISP20 = .0363, B2 = 
0.5, B3 = 0.5, NRADVL = 1 , A2 = 0.1, XF = 2.0, NCUT = 1 . RE = 0.8, WO = 
15.0, DIS = .05, NSC = 20. NSB = 8, NOS = 16, NSS = 2, NINPCS = 1 , BPO 
■ l/v/7, 0.0, etc., BML = 0.0. 

The numerical solution for the parallel-duct eigenfunctions and 
eigenvalues can be expected to become increasingly inaccurate as the 
spinning mode number increases. For m ^ 5 some caution should be used 
and for m ^ 10 numerical difficulties and inaccuracies are to be ex- 
pected. These high spinning modes will require an asymptotic expan- 
sion technique for accurate computation. 

5.8 Convergence of Results with Number of Modes 

It has been previously stated that one expects the results to be 
accurate if all cut-on modes are included in the calculation. In this 
section, a case for w = 18 is considered in which there are six cut-on 
modes in each direction. Calculations are made with fewer modes in 
order to check the convergence of the results as the number of modes is 
increased. The case discussed here is a straight duct with no mean flow 
and lined with a variable-admittance liner. Figure 11 shows the various 
approximations to the pressure profile at x = L/2 for the case when the 
lowest mode is input at x = 0. It is seen that convergence is slowest 
at the duct centerline, fastest at the wall, and that changes are quite 
pronounced across the entire duct for the first four modes. Addition of 
the fifth and sixth modes produces small changes (this would not be true 
if the input signal contained significant amounts of these modes). In 
contrast, the transmission coefficients converge more rapidly as shown 
in Table 1. With four modes, the transmission coefficients of all the 
modes that can be calculated are accurate. The coefficient of the 
lowest mode is the slowest to converge. Input data for calculations of 
these cases are W = 18.0, M = 0, IVK = 18.0, 17.59, 16.58, 14.85, 12.11, 



XMCO = 


7.26, -18.0, -17.59, -16.58, -14.85, -12.11, -7.26, ERR = 10"^, 

0.0, GAM = 1.4, RCF = 0.0, NTEMP = 3, DISP20 = .0363, B2 = B3 = 0.0, 
NRADVL = 0, A2 = 0.0, XF = 2,0, NCUT = 1 , RE = 0.0, BETAO = 0.0, BETAL = 
0.6 + 0.6i, NSC = 20, NSB = 2, NOS = 16. NSS = 2, NINPCS =1, BPO = 

1/v^, 0.0, etc., BML = 0.0. The value of NM changes from 2 to 12 as the 
number of modes is changed. 
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TABLE 1. Effect of increasing the number of modes on the 
convergence of the transmission coefficients. 


N 

tL,0 

'll 

jL.O 

'l2 

jL.O 

‘l3 

1 

-.323-.907i 



2 

-.251-.828i 

-.291-.349i 


3 

-.238-.825i 

-.309-.350i 

.080+.017i 

4 

-.236-.826i 

-.311-.348i 

.081+.012i 

5 

-.236-.827i 

-.311-.348i 

.081+.012i 

6 

-.236-.827i 

-.311-.348i 

.081+.012i 


-rUO 

‘l4 

-rL,0 

'l5 

tUO 

'l6 


1 

2 

3 

4 

5 

6 


-.012+.011i 

-.Oll+.Olli 

-.Oll+.Olli 


.001-.003i 

.001-.003i 


.0005+. 0001 i 



6. Summary 


An acoustic theory is developed to determine the sound transmission 
and attenuation through an infinite, hard-walled or lined circular duct 
carrying compressible, sheared, mean flows and having a variable cross 
section. The theory is applicable to large as well as small axial 
variations, as long as the mean flow does not separate. The technique 
is based on solving for the envelopes of the quasi-parallel acoustic 
modes that exist in the duct instead of solving for the actual wave, 
thereby reducing the computation time and the round-off error encountered 
in purely numerical techniques. The solution recovers the solution 
based on the method of multiple scales for slowly varying duct geometry. 

A computer program has been developed based on the wave-envelope 
analysis for general mean flows. The mean-flow model consists of a one- 
dimensional flow in an inviscid core and a quarter-sine profile in the 
boundary layer. Mean radial velocity effects can be included. Numeri- 
cal calculations performed for waves propagating in uniform ducts 
carrying fully-developed mean flows agree with the well-known results 
for uniform ducts. For non-uniform ducts, results are presented for the 
reflection and transmission coefficients as well as the acoustic pres- 
sure distributions for a number of conditions: both straight and vari- 
able area ducts with and without liners and mean flows from very low to 
high subsonic speeds are considered. The results for transmission and 
reflection coefficients are shown to possess symmetry characteristics in 
those cases for which it is expected. Comparisons with the results of a 
finite-element analysis for low-speed mean flows have shown reasonable 
agreement. Comparisons with one-dimensional results for high-speed mean 
flows have shown strong two-dimensional effects occurring near the duct 
throat. A number of test cases that demonstrate the flexibility of the 
program are included. Convergence of the transmission coefficients and 
of the acoustic pressure profiles with an increasing number of modes is 
illustrated. 

The only limitation of the wave envelope technique is that it is 
not suitable near cut-off, since the coefficient multiplying the term 
dA^/dx approaches zero. This problem is more apparent for a hard wall 
duct than for a soft wall duct, because k is exactly zero for a hard 


wall duct. Near cut off, the problem requires a turning-point analysis 

32 

using either the method of multiple scales or the Langer transformation . 
In addition, in a duct with very large axial gradients, such as occur in 
the throat region as M 1 , small axial steps may be required to assure 
that the program obtains an independent set of parallel-duct eigenfunc- 
tions to serve as basis functions. If two of the eigenvalues are the 
same an ill-conditioned matrix results and the program terminates. 



APPENDIX A 
Symbols and Notation 


Symbol 

32 Coefficient that specifies the minimum duct radius in Eq. 

(65) 

A Vector of mode amplitudes; 2N components 

A Amplitude of the n^^ right-running mode 

~ n iu 

A Amplitude of the n^" left- running mode 

n 

bi,b 2 Coefficients in the expression for the axial variation 
of the mean boundary-layer displacement thickness, 

Eq. (55) 

Vector of local amplitudes of the right- running modes; N 

4* h 

components; n^” component = A^exp{i/k|^dx} 

B" Vector of local amplitudes of the left-running modes; N 

components; n^^ component = A^exp{i/k^dx) 

c Speed of sound, c*/c* 

Cp Specific heat at constant pressure 

d Cavity depth of liner 

f An element of one of the coefficient matrices in the 

mn 

governing equation for the mode amplitudes, Eq. (50) 

F Coefficient matrix whose components are f^^ 

g An element of one of the coefficient matrices in the 

^mn 

governing equation for the mode amplitudes, Eq. (50) 

G Coefficient matrix whose components are g^^ 

i / -1 

th 

k^,k^ Complex propagation coefficients of the n right-running 
and left-running parallel-duct modes, respectively. 

L Duct length, l*/Rt 

m Circumferential mode number 

M Mach number of mean flow 

Mo Mach number of mean flow at x = 0 

N Number of parallel-duct modes propagating in each direction, 

p Pressure, p*/ptc*^ 

a a 

Pr Prandtl number, p*C*/k* 
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q An arbitrary physical variable 

f Position vector, r*/Rt 

r Radial co-ordinate, r*/Rt 

ri Recovery factor used in Eq. (54) 

R Duct radius, R*/Rl5f 

Rg Resistance of the liner facing sheet 

Re Reynolds number, ptc*Rt/y* 

ad W 

Ro Duct radius at x = 0 

R^’^ Reflection coefficient matrix at x = 0 for modes incident 
at X = 0; defined by Eq. (61) 

R^**" Reflection coefficient matrix at x = L for modes incident 
at X = L; defined by Eq. (61) 
t Time, t*c*/Ri5[ 

a 

T Temperature, T*/T* 

L 0 . ° 

T * Transmission coefficient matrix at x = L for modes inci- 

dent at X = 0; defined by Eq. (61) 

T^*^ Transmission coefficient matrix at x = 0 for modes inci- 
dent at X = L; defined by Eq. (61) 

Tg^ Adiabatic wall temperature; defined by Eq. (54b) 

TRi ,TR 2 jTRa ,TR 4 Transfer matrices defined by Eq. (60) 
u Axial velocity component, u*/c* 

a 

V Radial velocity component, v*/ct 

V Velocity vector, v*/c* 

a 

w Circumferential velocity component, w*/ct 

a 

X Axial co-ordinate, x*/R^ 

a Attenuation rate of the parallel-duct modes; imaginary 

part of k 

3 Liner admittance 

V Ratio of specific heats 

6 Mean boundary-layer thickness, 6*/Rt 

6i Mean boundary-! ayer displacement thickness, 5T/Rt 

6io Mean boundary-layer displacement thickness at x = 0 

n Variable defined in the adjoint homogeneous problem, Eqs. 

(44) - (47) 
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0 


Circumferential co-ordinate 

K Thermal conductivity, k*/k* 

W 

VI Viscosity coefficient, y*/u* 

w 

p density, p*/p* 

a 

T Viscous stress tensor, T*/p*c*^ 

<Pi »<p 2 t(p 3 94)u»4>5 9 (p 6 Weighting functions for the mass, axial 

momentum, radial momentum, circumferential momentum, 
energy, and state equations, respectively; related to the 
parallel-duct eigenfunctions through Eq. (48) and Eqs. 
(39)-{43) 

$ Viscous dissipation function, $*/p*c*^ 

a a 

\\) Parallel-duct eigenfunctions 

P T 

ip ,etc Parallel duct eigenfunctions for the acoustic pres- 
sure, temperature, etc. 

0 ) Circular frequency, 03*Rtf/ct = 27Tf*R^/ct where f* is the 

a a 

frequency in Hz. 

oa o) - kUo; Eq. (28) 

Subscri pts 

a Denotes a reference quantity; taken to be a mean-flow 

quantity at x = 0 and r = 0. 

0 Denotes a mean-flow quantity 

1 Denotes an acoustic quantity 

c Denotes a mean-flow quantity evaluated at the duct 

centerl ine 

mn Denotes a component of a matrix 

n Denotes either a quantity associated with a specific duct 

mode or an arbitrary component of a vector 

w Denotes a mean-flow quantity evaluated at the duct wall 

Superscripts 

( )^ Denotes a quantity associated with a right- running 
parallel-duct mode. 

( ) or ( ) Denotes a quantity associated with a left- running mode 

( )^ Denotes the transpose of a tensor 
* 

( ) . Denotes a dimensional quantity 



x=0 x=L 



Figure 1. Flow Configuration. 
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Figure 2. Comparison of wave-envelope and one-dimensional theories. 
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MAGNITUDE OF ACOUSTIC PRESSURE 


Figure 7. Variation of the acoustic pressure profiles with axial distance for a variable 
admittance and w = 9.12. 
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Figure lib. Effect of increasing the number of modes on the 
convergence of the acoustic pressure profiles at 
the throat. 
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